Spin dimers under staggered and random field in Cu2Fe2Ge40i3 



o 
o 

(N 
o 

o 

00 



I 



I 

o 
o 



> 
(N 
O 

G\ 
O 

o\ 

O 



T. Masuda/ K. Kakurai,^ and A. Zheludev^-'* 

'Department of Nanosystem Science, Yokohama City University, Yokohama, Kanagawa, 236-0027, Japai^ 
^Quantum Beam Science Division, JAEA, Tokai, Ibaraki 319-1195, Japan 
'^Institute for Solid State Physics, ETH, CH-8093 Ziirich, Switzerland 
^Laboratory for Neutron Scattering, ETH and Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland 

(Dated: October 28, 2009) 

We study 5 = 1/2 dimer excitation in a coupled chain and dimer compound Cu2Fe2Ge40i3 by 
inelastic neutron scattering technique. The Zeeman split of the dimer triplet by a staggered field 
is observed at low temperature. With the increase of temperature the effect of random field is 
detected by a drastic broadening of the triplet excitation. Basic dynamics of dimer in the staggered 
and random fields are experimentally identified in Cu2Fe2Ge40i3. 

PACS numbers: 75.10.Jm, 75.25.+Z, 75.50.Ee 



Excitations in quantum spin liquids can be viewed as 
strongly interacting bosonic quasiparticles. This circum- 
stance enables experimental studies of the physics of Bose 
liquids in prototypical quantum magnetic materials^'^''^. 
Such experiments are often possible under conditions 
that can not be realized in more conventional models, 
such as "^He^i^ and ultracold trapped ions^ii. One recent 
topic of interest is teh behavior of bosonic quasiparticles 
in the presence of disorder. Exotic new phases such as 
the Random Singlet stated, Bose and Mott glasses^ have 
been predicted for systems with quenched disorder. In 
real prototype materials one usually tries to create such 
disorder by chemical dopin g^^i^^ . In the present work 
we demonstrate an alternative approach: a random mag- 
netic field created by disordered (paramagnetic) ions. We 
show that such a random field acting on a simple dimer- 
based quantum spin liquid dramatically alters the exci- 
tation spectrum. 

Let us consider the effect of different types of mag- 
netic fields on an isolated S = 1/2 dimer, as shown in 
Fig. [TJ In a uniform field the excited triplet is spUt into 
three levels. Eventually, at high field, \S = 1, 5^ = 0) 
will cross the singlet ground state. In the presence of 
inter-dimer interactions, BEG of magnon will occur. If 
local fields applied to each dimer spin are antiparallel 
to each other (referred to as "staggered field" hereafter) 
the triplet is split into a singlet and a doublet. The sin- 
glet ground state becomes mixed with \S = 1,S^ = 0) 
and the total spin is no longer a good quantum number. 
Even in an infinitesimal staggered field the ground state 
becomes ploarized. Now, if the field direction is spatially 
randomized, each dimer will experience both a staggered 
and uniform component. The corresponding energy lev- 
els can be calculated numerically. The resulting density 
of state (DOS) for excitations in a set of N dimers is 
plotted in the right panel in Fig. [TJ The DOS lower and 
higher boundaries of the DOS spread coincide with the 
levels of \Sz ~ 1) and 15*2 = —1) in the uniform field. 

The quantum fcrrimagnct Cu2Fc2Ge40i3-^^ is a rare 
potential realization of this random field effect. The 
compound includes S — 1/2 Cu^+ dimers coupled to 
classical Fe'^+ chainsl2.. At low temperature the coop- 
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FIG. 1: Schematic diagrams of triplet excitations in S = 1/2 
dimers in different types of locally applied magnetic field. 



erative ordered state with classical spin and quantum 
spin is stabilized by a weak inter-subsystem coupling. In 
the adiabatic approximation, the quantum spins are ef- 
fectively under the internal field from the much slower 
fluctuating classical spins. In this compound, the stag- 
gered nature of the exchange field is due to the mag- 
netic structure. The staggered magnetization curves of 
dimers in Cu2Fe2Ge40i3^ were experimentally obtained 
by measuring the temperature dependence of sublattice 
moments in neutron diffraction. At high temperature, in 
the paramagnetic phase, the classical spins are thermally 
disordered and the effective field on the quantum spins is 
randomly oriented. Then the system can be considered 
as the ensemble of N dimers in a random quasi-static 
field. As shown in Fig. [1] the effect of this random field 
is to broaden the dimer excitations at T > Tn. 

In the previous inelastic neutron scattering study 
it was shown that the energy scales of excitations in 
the Fe chains and Gu dimers are well separatedi?-ii5'. 
The lower energy excitations up to 10 meV are Fe- 
based spin waves. Preliminary powder experimentsi^ 
and comparative studies in an isostructural compound 
Gu2Sc2Ge40i3i^ associated the dispersionless excita- 
tions at 24 meV with Gu-dimers. However the effect of a 
staggered and/or fluctuating field could not been identi- 
fied in powders samples. In the present paper we study 
the dimer excitations by single crystal inelastic neutron 
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scattering. By adopting a high resolution setup, we iden- 
tify the spHt peaks due to the staggered exchange field. 
Furthermore, we observe a drastic broadening of the peak 
profile at T > Tn that can be ascribed to randomly ori- 
ented field from thermally fluctuated Fc moments. 

High quality single crystals were grown by floating zone 
method. The crystal (monoclinic P2i/m) were found to 
be twinned, so that both microscopic domains share a* 
- b* plane. To avoid complications due to twinning, we 
restrict the measurements to the a* - b* plane. In the 
setups la and lb PG (002) were used for both monochro- 
mator and analyzer. The SoUer coUimations were 48' - 
60' - 60' - f20' and open - 80' - 80' - open for la and 
lb, respectively. In setup II, to achieve high energy res- 
olution, PG (004) for monochromator and PG (002) for 
analyzer with 30' - 20' - 40' - 120' were used. The setups 
la and II were performed on HBl spectrometer in HFIR, 
ORNL. The setup lb was performed on TASl spectrom- 
eter in JRR-3M, JAEA. In all setups final energy of the 
neutron was fixed at Ef = 14.7meV and PG filter was 
installed after the sample to eliminate higher order con- 
tamination. A closed cycle He refrigerator was used to 
achieve low temperatures. 

In a series of energy scans in a wide range of {h k 0) 
space shown in Fig[2{a) two dispersionless peaks are read- 
ily identified: a pronounced one at fto; 24 meV and a 
weaker feature at fio; 31 meV. The experiments were 
performed in setups la and lb. The former is consis- 
tent with the Cu-centered magnetic excitation in previ- 
ous studiesi^ii^. Constant energy scan at fto; = 24 meV 
and its temperature dependence are shown in FigllKb). 
The observed sinusoidal intensity modulation is charac- 
teristics of dimer excitations and is observed in a wide 
temperature range. In fig. [UJc) the temperature depen- 
dence of the peak intensity is shown. The intensity at 
q = (0 2.5 0) was measured at each temperature and then 
subtracted as background. The decrease of the intensity 
at high temperature is common behavior for magnetic 
excitations in local spin clusters. The smaller peak at 
hiUj ~ 31 meV was identified as a Fe-centered excitation, 
as will be discussed below. 

To obtain a more detailed profile, we performed an 
energy scan using setup II at T = 2.0 K. As shown in 
Fig. [3] it is revealed that the primary peak at fiw = 24 
meV actually has a shoulder structure. The main peak 
is located at 23.5 meV, and a smaller bump is centered 
around 25.0 meV. This splitting is attributed to the stag- 
gered exchange field from the adjacent Fe moments. The 
main peak corresponds to the excitation doublet and the 
small one to the singlet. 

Energy scans collected at several temperatures are 
shown Fig. [DJa). The small peak at /iw ~ 31 meV in 
Fig. [5] is temperature independent and has been sub- 
tracted from the data. Well-defined peaks are observed 
at all temperatures. While at low temperature the peak 
profile is sharp and the width is within resolution limit, 
at T > T/v the peak becomes drastically broadened. 
This qualitative behavior is consistent with the effect a 
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FIG. 2: Inelastic neutron scattering using experimental setup 
la and lb. (a) Typical energy scans at (h k Q). Dispersionless 
excitations are observed at huj = 24 and 31 meV. Two peaks 
are separately fit by Gaussians (dotted curves), (b) h scans at 
hjjj = 24 meV at various temperatures. Sinusoidal intensity 
modulations are fitted to the dimer structure factor calculated 
for zero field, plus a constant background (solid curves), (c) 
Temperature dependence of the peak intensity at q = (/i 2.5 0) 
and hw = 24 meV. 
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FIG. 3: Energy scan collected using the high resolution setup 
II at T = 2.0 K. The shoulder structure is reproduced by the 
doublet (shaded with white background) and singlet (shaded 
with gray background) dimer excitations, split by a staggered 
exchange field. 



random exchange field should have on the dimer exci- 
tation triplet. The data were analyzed using Gaussian 
fits. The estimated peak positions, widths, and the inte- 
grated intensities are plotted as functions of temperature 
in Figs, mb)-(d). With increasing temperature the peak 
energy decreases at T Tn, and stays constant beyond. 
The peak width drastically increases at T ~ Tn, but also 
remains constant at higher temperature. The integrated 
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FIG. 4: (a) Temperature dependence of the peak profile at 
q — {2 3 0). Experimental resolution is indicated by the gray 
area. Small peaks due to Fe-centered excitation at 31 meV 
are separately fitted and subtracted. Profiles at T > 54 K 
are reproduced by dimers in a randomly oriented field model 
(solid curves) . The temperature dependence of peak positions 
(b), widths (c) and integrated intensities (d), as estimated 
from Gaussin fits. 



intensity decreases by 10 ^ 20%. It is noted that in the 
previous powder experiment the peak cannot be distin- 
guished at T > 41 Ki^. This is because the powder inte- 
gration in wide q space collects phonon excitations and 
accidental suprious peaks, masking magnetic excitations 
at higher temperatures. 

For T < Tn we will consider the following effective 
Hamiltonian: 

H = JcuSi ■ S2 + .g/iB'S'i • hi + gnBS2 ■ ^2, (1) 

where hi — {0 h) and /i2 = (0 — h). Here the z- 
axis is chosen along the ordered Cu moment. The ground 
state energy Eq decreases with the field, Eq — —Jcu/4: — 
\J[2g^'QhY + and the excitation triplet splits into 

a singlet and doublet. The corresponding energy levels 
are giben by 



(2) 



Ad = Jcu/2 + ^{2gmhY + JlJ2, (3) 
and plotted vs. h in Fig.[TJ For gji-eh A, the neutron 



cross section is approximately given by: 

-^^A.(,.„)^^sin^(,.d)(/(,))^P(T) 
x{A(/i)(l + cos2 e)6{hLj - Ad)+B{h) sin^ 0d{huj - A,)}.(4) 

The doublet and singlet terms correspond to transverse 
and longitudinal spin fluctuation, respectively. A{h) and 
B{h) are h dependent parameters with A{h) < 1, B{h) < 
1, and A{0) = B{0) = 1. Since the staggered field sta- 
bilizes the polarized spin configuration and suppresses 
longitudinal fluctuation, B{h) decrease with h. Mean- 
while A{h) is almost constant in the low fleld. P{T) 
is a temperature factor, P{T) — 1/{1 + 2exp(— /3Ad) -I- 
exp(— /3As}. q is the scattering vector, d is the spin sepa- 
ration in each dimer, and 9 is the angle between q and the 
moment of Cu. We used mcu = (-0.227, 0.035, -0.301) 
/iB^^ to calculate 9. Two types of domains, namely anti- 
fcrromagntic and crystallographic ones due to twinning, 
are considered. 

The peak profile in Fig.[3]is reasonably well reproduced 
by the cross section convoluted by experimental resolu- 
tion function with Ag — 25.0 meV and Ad = 23.5 meV. 
From eqs.dl]) and (O, Jcu = 22.0 meV and h = 51 T are 
obtained. Let us check the consistency of h with the pre- 
vious study From the staggered magnetization curve 
by neutron diffraction Jcu-Fc/>/cu = 0.105 was obtained. 
Here Jcu-Fc is the interaction between Cu and Fe spins. 
Using the molecular field relation h = TOFc^Cu-Fe/(<?AiB)^ 
and previously obtained parameters, ft, ~ 40 T is esti- 
mated. Thus statically estimated value is consistent with 
that obtained in the present dynamic measurement. 

The energy splitting between the singlet and doublet 
states is about 1.5 meV. This value is small compared 
energy resolution in the typical experimental setup. In 
setup la and lb at T < Tn, therefore, staggered field ef- 
fect is smeared and two terms in eq.(|4]) are integrated. 
Then the cross section is approximately equivalent to 
that at /i = 0. Indeed, the constant energy scan at T 
= 3.3 K in Fig.jSJb) is reasonably fitted by dimers cross 
section in zero field shown by the thick curve. 

At r > Tn effective field on the Cu dimers is ran- 
domly oriented. We consider an ensemble of N dimers in 
random field. The randomly oriented field h is assumed 
to have a constant magnitude in eq.(IT]). The resulting 
DOS of the excited states is then calculated numerically. 
The neutron cross section is assumed to be approximately 
proportional to the DOS, 
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= {^rof'^sm^q.d)f{qfP,,,,^{T)D{nu) (5) 



with / D{e)de = 3N and Prand(T) = N/{N + 

J D{e)e~'^de). The data collected at T > 54 K are wefl 
reproduced by this cross section convoluted by experi- 
mental resolution function, as indicated by solid curves in 
Fig. m^a). The obtained fit parameters are Jcu — 22.3(4) 
meV and h = 41.(8) T. The values are reasonably consis- 
tent with those obtained at T < Tn. The q dependence 
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of the cross section is the same as for zero field, and is 
given by dimer structure factor sin^(q.rf). Indeed, the q 
scans at 80 K and 300 K in Fig. [SJb) are reproduced by 
this model. The temperature dependence in Fig. is 
well accounted for by the temperature factor Prand(T)- 

We shall now discuss the small decrease of the inten- 
sity at T < Tn in Fig. [IJd). At T > Tn dimer spins 
are fluctuated equally in all directions and the dynami- 
cal spin correlation is fully detected by neutron. In the 
ordered state the polarized magnetic ordering suppresses 
the longitudinal fluctuation of Cu spins. To estimate 
the reduction of the longitudinal excitation we calculate 
B{h = 51T) = 0.77. The reduction of is about 20% 
that is consistent with the experiment. This means that 
51 T is rather modest compared with the intradimer in- 
teraction Jcu = 22 meV. If the effective field was large 
and the moment were fully polarized, the suppression 
would be more drastic. Such a situation is in fact re- 
alized in Haldane spin chains coupled to rare earth mo- 
ment in Pr2BaNi05 with fully saturated Ni^+ moment 
at T < T^r^. The Haldane-gap mode lost half of its 
intensity at T < Tn and it was ascribed to the total 
suppression of longitudinal mode. 

Finally we will mention the temperature independent 
small peak at ftw - 31 meV in Fig.EIJa). If the Fe 5" = 5/2 
chains were perfectly isolated from the Cu subsystem, 
the Fe excitation spectrum would be dominated by one- 



magnon excitation at huj < 5Jfc- However, a recent the- 
ory predicts that the introduction of Cu dimer enhances 
the multi-magnon excitation of Fe spins at huj = IOJfc, 
15Jfo, 20Jfc and 25Jfc. According to the Bond operator 
method^'*' the excitation at huj = 20Jfc is the partic- 
ularly enhanced^S. Since Jpc = 1.6 meV^^, the observed 
small peak at ftw = 31 meV could be ascribed to the Fe 
centered longitudinal excitation. Further details will be 
published somewhere else. 

To conclude, we have experimentally investigated the 
dynamics of S* = 1/2 dimers in staggered and random 
fields in Cu2Fe2Ge40i3. The staggered field is realized at 
T < Tn and produces a splitting of the excitation triplet. 
At T > Tn a random exchange field produces a drastic 
broadening of these modes. In teh future, polarized neu- 
tron experiments may be useful to separate the longitu- 
dinal and transverse excitations. Recently Cu2CdB20g^ 
and Cu3Mo20g^^ identified as new realizationsof the cou- 
pled dimers and chains models. Particularly in the latter 
compound, the dimer energy is close to that of the chains, 
and more complex physics is expected. 
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